PROCEEDINGS 

™ SCIENCE 



Perturbative calculations for the HISQ action: the 
gluon action at 0(Nfasa^) 



A. Hart^ G.M. von Hippel*^, R.R. Morgan^ 

" SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3JZ, U.K. 
* Deutsches Elektronen-Synchrotron DESY, 15738 Zeuthen, Germany 
" DAMTP CMS, University of Cambridge, Cambridge CBS OWA, U.K. 

Email: georg . von . hippelSdesy . de 



We present a new (and general) algorithm for deriving lattice Feynman rules which is capable 
of handling actions as complex as the Highly Improved Staggered Quark (HISQ) action. This 
enables us to perform a perturbative calculation of the influence of dynamical HISQ fermions 
on the perturbative improvement of the gluonic action in the same way as we have previously 
done for asqtad fermions. We find the fermionic contributions to the radiative corrections in the 
Luscher-Weisz gauge action to be somewhat larger for HISQ fermions than for asqtad. 
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1. Introduction 

Continuing rapid advances in parallel computing, along with theoretical progress in the formu- 
lation of lattice field theories with fermions, have led to lattice QCD simulations with dynamical 
light quarks becoming the norm rather than the exception. 

The Fermilab Lattice, MILC and HPQCD collaborations have an ambitious program which 
to date has made several high-precision predictions from unquenched lattice QCD simulations [1]. 
This body of work is based on the Symanzik-improved staggered-quark formalism, specifically the 
use of the asqtad [2] action. More recently, the Highly Improved Staggered Quark (HISQ) action 
has been used to further suppress taste-changing interactions and to allow the use of heavier quarks 
at the same lattice spacing by removing tree-level 0((ma)^) artifacts from the quark action [4]. In 
order to consistently use the HISQ action for the sea quarks as well [5], the calculation of HISQ 
quark loops on the Symanzik-improvement of the gluon action is also needed. Having previously 
carried out that calculation for the asqtad action [6], we update our calculation here to apply to the 
case of dynamical HISQ fermions. 

2. Perturbation Theory for the HISQ action 

The HISQ action is defined by an iterated smearing procedure with reunitarisation: 

U"''Q = (Fasq' O Pj,(3) O Fpatv) [U] (2. 1) 

where U = exp(gA) is the unsmeared gauge field, Pu{3) denotes the polar projection onto U{3) 
(as used in simulations, and not SU{3)), and the Fat7 and modified asq smearings are defined in 
[4]. Straightforward application of the methods from [7] to this action is unfeasible, since the 
memory requirements for expanding the action directly into monomials quickly become excessive. 
We therefore take advantage of the two-level structure inherent in the definition of the action and 
split the derivation and application of the Feynman rules into two steps. 

In the first step, the Feynman rules for the outer layer (the modified asqtad action) are derived 
in the same way as previously. We use our HiPPy python code [7] to expand the asq' action in 
terms of the Fat7R smeared link 

[/P^™(x) = (Pj;(3) oF^,a)[U,] = A(-+2A) (2.2) 

with a Lie-algebra-valued field , giving the usual monomials 

= fr ¥{Xr,i)B^,, {Vr;n) ■ ■■B^,XVr-i,r)'rrMyrj) (2-3) 

To derive the full HISQ Feynman rules, we also need to know the expansion of in terms of the 
original gauge potential . To obtain this, we write the Fat7-smeared link as ' Fpatv [U] = M = HV, 
where = H and V G U{3). We can now use our HiPPy expansion routines [7] to obtain an 
expansion 

M = c[l + an *A^+a^v * {A^ *Av) + . . .] (2.4) 
'in the following, we will suppress Lorentz and lattice site indices. 
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where, e.g. a^v*(A^*^v) ='Lx.y<^iiv{^,y)Afi{x+2fi-)Av{y + 2V). Then unitarity of V implies that 
7? = MM^ = and hence V = R^^I^M using the expansion 

R-"^ = {\ + {R-\))-"^ = \-\[R-\) + \{R-\Y + ... (2.5) 
Rearranging the result as V = exp(B), i.e. 

B = \og{y) = {v-\)-^-{y-\f + ... (2.6) 

finally yields the desired expansion of B. 

Given this, we can now numerically reconstruct the HISQ Feynman rules for any given set of 
momenta from eqn. (2.3) by a convolution of the asq' Feynman rules of eqn. (2.3) with the expan- 
sion of in terms of A^, summing up all the different ways in which the gluons going into the 
vertex could have come from the fields appearing in eqn. (2.3). Compared to a simple-minded 
expansion of the HISQ action, this not only save enough memory to enable the derivation to be per- 
formed in practice, but also leads to a considerable speed-up in many cases. In particular, we can 
take advantage of the (anti-)symmetries that the expansion of in terms of possesses, allowing 
us to reduce the number of contributions we need to take into account when evaluating Feynman 
diagrams. In the calculation of the three-gluon vertex for the "octopus" diagram (a fermion tadpole 
with thi^ee gluon legs) entering the three-point function, we are able to omit the contribution from 
the expansion of a single into three gluons on symmetry grounds. 

3. On-shell improvement 

The Liischer-Weisz action is given by [8] 

■5 = 11(1-8(^1+^2)) £ (1-Pmv>+2ci £ (l-/?;,v> + ^C2 £ (l-7;,vp>[, (3.1) 

where P, R and T are the plaquette, rectangle and "twisted" parallelogram loops, respectively. The 
coefficients c\ and C2 need to be determined in order to eliminate the 0'{a^) lattice artifacts. 
Given two independent quantities Q\ and Q2 with expansions 

Qi = Qi + Wiinaf + dijCjinaf + ^ ((Ata)^) , (3.2) 

in powers of (/ia), where jJ. is some energy scale, we obtain the i^(a^) matching condition 

dijCj = -Wi . (3.3) 

Since this equation is linear, both sides can be decomposed into a gluonic and a fermionic part; the 
known gluonic part [9, 10] being independent of the fermion action, we will here focus only on the 
fermionic part. 

At tree-level, there are no fermion loops to consider, and hence the tree-level coefficients re- 
main unchanged compared to the quenched case [9]. To compute to one-loop fermionic corrections 
to the gluon action, we will follow the same procedure as in the case of the asqtad action [6]. 
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4. Twisted boundary conditions 

We work on a four-dimensional Euclidean lattice of length La in the x and y directions and 
lengths L^a, Lta in the z and t directions, respectively, where a is the lattice spacing and L,L^,Lt 
are even integers. In the following, we will employ twisted boundary conditions in much the same 
way as in [9, 10]. The twisted boundary conditions we use for gluons and quarks are applied to the 
{x,y) directions and are given by (v = x,y) 

u^{x+Lv)=ayU^{x)a-^ , ^'{x+Lv)=ay^>{x)a-^ , (4.i) 

where the quark field ^sc{x) becomes a matrix in smell-colour space [11] by the introduction of a 
new SU(N) quantum number "smell" in addition to the quark colour. We apply periodic boundary 
conditions in the {z,t) directions. 

These boundary conditions lead to a change in the Fourier expansion of the fields: in the 
twisted {x,y) directions the momentum sums are now over 

NL NL , ^ 

py=mny, - — <ny < —, v = {x,y) , (4.2) 

where m= j^. The modes with {rix = ny = mod A'^) are omitted from the sum in the case of the 
gluons. The momentum sums for quark loops need to be divided by A'^ to remove the redundant 
smell factor. 

The twisted theory can be viewed as a two-dimensional Kaluza-Klein theory in the (z, t) plane. 
Denoting n = {nx,ny), the stable particles in the (z,f) continuum limit of this effective theory are 
called the A mesons (n = ( 1 , 0) or n = (0, 1 )) with mass m and the B mesons (n = ( 1 , 1 )) with mass 
V2m [10]. 



5. Small-mass expansions 

Although we ultimately wish to extrapolate to the chiral limit, we cannot set m^ja = straight 
away, since the correct chiral limit is triqa — > 0, ma 0, niq/m > C, where C is a constant de- 
termined by the requirement that a Wick rotation can be performed without encountering a pinch 
singularity. 

Therefore, we first expand some observable quantity Q in powers of ma at fixed m^a: 

Q{ma,mqa) = a^^\mqa) + a^.^\mqa){ma)^ + & ((ma)^, (ma)^log(m<3)) (5.1) 

where the coefficients in the expansion are all functions of m^a. There is no term at ((ma)^ log(ma)) 
since the gluon action is improved at tree-level to 0{a^) [10]. Then, we expand the coefficients 
a^^\mqa) in power of mqa. 
For a^^\mqa) we have 

a^^\mqa) = b^^^\og{mqa) + af^ . (5.2) 

Since we expect a well-defined continuum limit, a^^\mqa) cannot contain any negative powers of 
mqa, but may contain logarithms; b^^^ is the anomalous dimension associated with Q, and can be 
determined by a continuum calculation. 
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□ data 
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Figure 1: A plot of the fermionic contributions to the one-loop A meson self-energy mj^^'/w against (ma)^. 
The vanishing of m'^^ /m in the infinite- volume limit can be seen clearly. 

For a'^\mqa) we find 

af\mqa) = +afi + (4§ +bfl log(m^a)) {mquf + ^ (Ka)^) . 

After multiplication by [maf-, the (niqa)^^ contribution gives rise to a continuum contribution to Q, 
and ^2^2 calculable in continuum perturbation theory. There can be no term in {niqa)^^ \og{mqa) 
since this would be a volume-dependent further contribution to the anomalous dimension of Q, and 
there can be no term in log(mga) since the action is tree-level 0{a^) improved. 

In the chiral limit niq 0, the term w; that appears on the right-hand side of Eqn. (3.3) is o'^q. 

6. Twisted spectral quantities 

The simplest spectral quantity that can be chosen within the framework of the twisted boundary 
conditions outlined above is the (renormalised) mass of the A meson. The one-loop correction the 
the A meson mass is given by 



'a --2o(k)— ^ 
2m\' 



(6.1) 



k=[im^°\0,mfi) 



where Zo(k) = 1 + ^ {{maY^ is the residue of the pole of the tree-level gluon propagator at spatial 
momentum k, and m^^^ is defined so that the momentum k is on-shell. 

From gauge invaiiance we find = and aQ"'''^'(m^a) = 0. The (^as{ma)^) contribu- 

tion from improvement of the action is given by [10] 

(1) 

Aimp^ = -(c'/'-4'^)W + ^(M') • (6-2) 
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Figure 2: Plots of Oq against niqa (left) and of against niqa (right) with the fits shown for compari- 



son. 



The next simplest independent spectral quantity is the scattering amplitude for A mesons at B 
meson threshold, which can be described by an effective AAB meson coupling constant A [12]: 



X=goVZ{^)Z{p)Z{(i)ejr''^''{k,p,q) 



(6.3) 



with a twist factor of ^Tr([ri, PpjE^) factored out from from both sides, and the momenta and 
polarisations of the incoming particles are (where r > is defined such that E{q) = 0) 



k={iE(k),k) p = {-iE{p),p) q = iO,q) 



e = (0,l,-l,0) 



k=(0,m,/r) p = (m,0,/r) q = (— m, — m, — 2/r) 
We expand Eqn. (6.3) perturbatively to one-loop order and find (up to ff{{nia)^) corrections) 



(6.4) 



m 



1 



1 ,\ri'^ 
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ko=iE{k) 
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12 " J dq} V 



e'e^nipiq) 



qo=0 



(6.5) 



The derivatives of Feynman diagrams are computed using automatic differentiation [13]. Contin- 
uum calculations of the anomalous dimension and infrared divergence give 

The improvement contribution to A is [10] 



^2,-2 



120;r2 



(6.6) 



.(1) 



-»imp 



m 



(6.7) 



7. Results 



To extract the improvement coefficients from our diagrammatic calculations, we compute the 
diagrams for a number of different values of both L and with Nf = I, N = 3. At each value 
of niq, we then perform a fit in ma of the form given in Eqn. (5.1) to extract the coefficients 
ai^'^\mqa), n = 0,2. Our fits confirm that al!^"^'^\mqa) = 0. 
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Performing a fit of tlie form (5.2) and (5.3), respectively, on tliese coefficients, we are able 
to extract the analytically -known coefficients with high accuracy, along with the required (ma)^ 
contributions. 

Solving equation (3.3) for c|'' given the fitted values for ajQ. our results can be summarised 

as 



where the quenched (Nf = 0) results are taken from [10]. The shift from the unquenched values 
is surprisingly large, even compared to the coefficients for asqtad fermions [6]. At first sight, this 
may seem like a surprise, since HISQ is supposed to be the more highly-improved action. However, 
HISQ is designed to suppress taste-changing interactions (low momentum quark/high momentum 
gluon couplings), but these coefficients come from high momentum quark/low momentum gluon 
couplings, for whose suppression the HISQ action is not tuned. 
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